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ABSTRACT 

We model the optical to far- infrared SEDs of a sample of six type-1 and 
six type-2 quasars selected in the mid-infrared. The objects in our sample are 
matched in mid-IR luminosity and selected based on their Spitzer IRAC colors. 
We obtained new targeted Spitzer IRS and MIPS observations and used archival 
photometry to examine the optical to far-IR SEDs. We investigate whether the 
observed differences between samples are consistent with orientation-based uni- 
fication schemes. The type-1 objects show significant emission at 3 /xm. They 
do not show strong PAH emission and have less far-IR emission on average when 
compared to the type-2 objects. The SEDs of the type-2 objects show a wide 
assortment of silicate features, ranging from weak emission to deep silicate ab- 
sorption. Some also show strong PAH features. In comparison, silicate is only 
seen in emission in the type-1 objects. This is consistent with some of the type-2s 
being reddened by a foreground screen of cooler dust, perhaps in the host galaxy 
itself. We investigate the AGN contribution to the far-IR emission and find it to 
be significant. We also estimate the star formation rate for each of the objects 
by integrating the modeled far-IR flux and compare this with the SFR found 
from PAH emission. We find the type-2 quasars have a higher average SFR 
than the type-1 quasars based on both methods, though this could be due to 
differences in bolometric luminosities of the objects. While we find pronounced 
differences between the two types of objects, none of them are inconsistent with 
orientation-based unification shemes. 
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1. Introduction 

The study of various types of AGN can provide insight into their origin and evolution. 
AGN are classified according to their optical emission lines as type-1 or type-2. Type-1 AGN 
have an "unobscured" view of the hot accretion disk and broad emission lines originating 
from fast moving clouds that surround a supermassive black hole. Type-2 objects are known 
as "obscured" AGN, because no broad lines or continuum emission from the accretion disk 
are observed, yet the spectra still show narrow emission lines similar to those of the type-1 
objects. 

The interpretation of the differences between type-1 and type-2 quasars (the most lumi- 
nous AGN) remains an open issue. The orientation-based unification scheme proposes the 
two types of objects are intrinsically the same, but appear different observationally due to 
the presence of an obscuring torus of dust close to the central power source (see Antonucci 
1993). Type-1 objects are seen from an orientation which allows an unobscured view of the 
central source. This allows the broad line region of the quasar to be observed. However, 
type-2 objects are seen through the dusty torus material, which hides the central contin- 
uum source and broad-line region. Evidence that this model is at least partly correct comes 
from observations of polarized broad-line emission in galaxies classified as type-2 quasars 
(e.g., Zakamska et al. 2005). But this model does not provide insight into the origin of 
quasar activity. 

An alternative model propose s that quasar activit y is triggered by galactic mergers and 



evolves through an obscured stage fISanders et al.lll988l ). As the galaxies merge, gas and dust 



are stirred up, triggering quasar activity and star formation. Initially the quasar would be 
obscured by the dust, but over time will become unobscured as dust and gas are blown away 
from the nuclear region by quasar and/or supernovae driven winds and radiation pressure. 
Indeed, several host galaxy studies have been made which show that luminous AGN are 
hosted by massive, mainly spheroid-dominated galaxies, with a significant fraction, but not 
all, showing signs of a relatively recent merger and associated star formation (e.g., Dunlop 
et al. 2003; Canalizo et al. 2007; Bennert et al. 2008). Thus the evolution model predicts 
that the host galaxies of type-2 quasars should have more recent star forming activity than 
those of type-Is. Although the discovery of polarized broad lines is a powerful argument in 
favor of the orientation-based model, it does not rule out an evolutionary link. 

The obscuring dust in a host galaxy would be apparent in its observed SED. Dust close 
to the AGN would be heated by the central power source and would thermally re-radiate 
in the mid-IR (100 - lOOOK), dominating the rest-frame SED from ~2 - 40 /xm. This high 
temperature dust emission makes the mid-IR SEDs of both optically-obscured and optically- 
unobscured AGN quite distinct from those of either normal or starburst galaxies, whose dust 
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temperatures are always < lOOK. By modeling the various components of the SEDs, we can 
compare the relative amounts of "cool", "warm", and "hot" dust emission present in the 
host galaxies. This will give us some idea of what mechanism dominates the heating of the 
dust (star formation or AGN), and how much of a role orientation might play in the observed 
SEDs. 

In this paper, we present modeled fits to the SEDs of a sample of type-1 and type-2 



quasars of similar redshift and luminosity. iLacy et al.l (120071 ) presented preliminary fits to 



the type-2 quasars in this paper. We also measure the objects' far-IR luminosities, PAH 
luminosities, and estimate their star formation rates. We then compare the results with 
starburst galaxies. We describe the sample selection in §2, and the observations and re- 
ductions in §3. In §4 we present the analysis and results including our modeling procedure 
and resultant spectra, a color-color diagram, and star formation rate (SFR) calculations. 
Finally, we summarize our results in §5. Throughout this paper we adopt the cosmological 
parameters of Ho=71 km s~^ Mpc^^, ^7^=0. 27 and Q\=0.73. 



2. Sample Selection 



Quasar surveys are now highly effective at finding quasars on the basis of their optical 
colors, provided they are not reddened or only lightly reddened by dust. With surveys such 
as the Sloan Digital Sky Survey (SDSS) (Schneider et al. 2003), 2dF (Groom et al. 2004), 
and Gombo-17 (Wolf et al. 2003), a census of such objects across the full range of quasar 
luminosities and redshifts out to z ~ 6 is largely complete. Only recently have radio-quiet 
type-2 quasars been found in substantial numbers. However, there is still no consensus on 
how the dust-obscured quasars are related to the unobscured population. In particular, it is 
difficult to select matched samples of type-1 and type-2 objects with which to perform such 
studies. By using Spitzer data from the extragalactic First Look Survey (XFLS), we have 
been able to construct samples of type-1 and type-2 quasars matched in mid-IR luminosity 
and redshift. 



Lacy et al.l (l2004j ) showed that quasars occupy a distinct "sequence" of a Spitzer IRAG 
color-color plot (their Fig. 1). Follow-up spectroscopy of a sample of mid-IR selected AGN 
and quasar candidates showed that about 2/3 of AGN selected in this manner may be 
sufficientl y obscured by d ust in the optical to be missing from optically-selected quasar 
samples (ILacy et al.l 120071 ). We selected six type-2 quasars an d six t ype-1 quasars (see 
Table [1]) from the 1 mJy 8 /im flux- limited sample of iLacy et al.l (120041 ). with similar mid- 
IR colors and 5.8 /im luminosities. These objects are approximately matched in mid-IR 
luminosity. Our modeling (see §4.1!) shows they have mid-IR power-law luminosities within 
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a factor of five from each other. Furthermore, the observed frame 24 /xm luminosities of the 
type-1 and type-2 quasars are not statistically different from each other. It is important 
that the sample be matched in luminosity, since the quasar luminosity will affect the rate at 
which dust and gas is removed from the nuclear region. 



Observations and Data Reduction 



We supplement dat a from the S'pztzer Extragalactic First Look Survey (xFLS) (iFadda et al 



20061 : iFrayer et al.1120061 ) with targeted observations in the 70 /im and 160 /im using the Multi- 
band Imaging Photometer for Spitzer (MIPS). Each of our targets was detected at the 24 
lim wavelength in the Spitzer xFLS. Five of our targets were detected in the 70 fim xFLS 
mosaic, and we measured the 70 /xm fluxes for two other targets directly from the xFLS mo- 
saic image. The remaining five targets of our sample were not detected in the 70 /im xFLS 
mosaic and so required targeted observations. We observed each of these targets for 500 
seconds. All of the targets were observed with the MIPS 160 fj.m band for 2000 seconds 
each. 

The photometric data were reduced through the standard Spitzer MIPS pipeline. We 
used the post basic calibrated (post-BCD) filtered mosaic images to measure the flux of each 
target. The pixel scales of the post-BCD MIPS images are 4 arcsec pixel"^ at 70 /xm and 8 
arcsec pixel"^ at 160 /im . 

We also observed each target using the short-low and long-low resolution modes with 
the Spitzer infrared spectrograph (IRS). The spectra were reduced through the standard 
Spitzer pipeline. We manually masked the bad pixels in the background subtracted spectra 
(difference between nodded exposure positions) using the IDL program IRSCLEA]N(^. We 
then extracted the spectra using the standard point source extraction in the Spitzer SPICeI^ 
software tool, and combined the individual spectra to create a continuous mid-IR spectrum 
for each target. Lastly we included photometry measured from our MIPS data (Tabled]) and 
data available for our targets from SDSS, 2MASS and Spitzer IRAC observations (Table [2]) 
to produce an observed frame SED that ranges from optical through the infrared. 



^http:/ /ssc. spitzer. caltech.edu/archanaly/contributed/irsclean 
^http: / / ssc.spitzer.caltech.edu / postbcd/spice.html 
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3.1. Photometry 

We used the DAOPHOT package in IRAF to measure the 70 fim and 160 fim fluxes 
of our targets. For the 70 /xm measurements we examined each image individually and 
determined the best aperture and sky annulus that maximized the signal-to-noise ratio. We 
then applied aperture corrections from the Spitzer Space Telescope MIPS Data HandboolJ^, 
and compared the measured flux to a 3a limit obtained by measuring the noise in the image. 

Each of our 12 targets was observed in the 160 fim MIPS band, though many of our 
targets were not detected at this wavelength. We chose an example image with a strong 
detection and determined the aperture and sky annulus sizes that maximized the signal- 
to-noise ratio. We used those parameters while measuring the other images with bright 
detections. The same procedure was used for the images with targets that were not detected. 
We then compared the measured fluxes of each of the targets to the 3a limit determined 
from the noise in the image. Five of the 12 targets have fluxes less than the 2a noise limit, 
and we count them as non-detections. Five other targets have fluxes >2a, and only one has 
160 /im flux >3a. The remaining target is discussed following. 

The target SDSS 1711-1-5855 required special attention, because there was a bright 
object near the target. We flt a PSF using several other bright objects in the image, then 
subtracted sources from the image leaving only our target. We then measured the flux of 
SDSS 171H-5855 from the subtracted image. This procedure was successful at 70 /im, but 
our attempts at removing the extra source at 160 /xm resulted in an image with signiflcant 
residuals. Therefore, we only report an upper limit measured from the combined flux of the 
objects at 160 fim. This upper limit was not used in our modeling procedure and analysis. 

Table [1] shows the results of our MIPS flux measurements combined with data from the 
xFLS. 



3.2. Spectra 

We present the Spitzer IRS spectra for our sample in Fig. [H The spectra show distinct 
differences. The spectra of the type-1 quasars tend to have fewer features and are more 
homogeneous than those of the type-2s. The type-1 objects show signiflcant continuum 
emission near 3 //m, little PAH emission, and no signiflcant silicate absorption. The type-2 
spectra are more heterogenous, showing a variety of silicate absorption depths and PAH 



•^http:/ /ssc. sptizer.caltech.edu/mips/apercorr/ 
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strengths, consistent with previous observations of AGN (e.g., Hao et ah 2007; Spoon et 
ah 2007; Zakamska et ah 2008) as well as starburst galaxies, ULIRGs and sub-mm galaxies 
(e.g., Brandl et al. 2006; Sajina et al. 2007; Lutz et al. 2008; Menendez-Delmestre et al. 2009). 
The type-2 quasars have on average double the PAH luminosity of the type-1 quasars (see 
Table [3]). 



Spoon et al.l (120071 ) classify spectra based on their 6.2 /im PAH equivalent width and 



9.7 /im silicate absorption strength. Our type-2 spectra are similar to their IB, IC and 
2B classifications. These classifications are dominated by Seyferts, starbursts, and ULIRGs, 
respectively. Our type-1 quasars have spectra most similar to the lA classification of galaxies 
by Spoon et al. 

The variation of silicate abs orption strength s in th e type-2 spectra of our sample is 
also consistent with that seen by IZakamska et al.l (120081 ). who study a sample of optically 
selected type-2 quasars from the SDSS. Zakamska et al. find silicate absorption similar to 
some of our sample, including two objects that have more significant absorption than our 
sample (see Lacy et al. 2007). Four of their objects appear to have no silicate absorption, 
including one that may show silicate emission, consistent with our mid-IR selected type-1 
quasars. They al so address the correlation of silicate absorption with axial ratio found by 
Lacy et al.l (120071 ). Lacy et al. interpreted this correlation to indicat e that the silicate dust 
resides in the host galaxy, and not the nuclear region. The objects of I Zakamska et al.l (120081 ) 
show no correlation, and they propose this could be a difference in samples due to selection 
biases based on luminosity. 



4. Analysis and Results 
4.1. Modeling 



We modeled the SEDs of our sample from the optical out to the far-IR. We used avail- 



able SDSS and 2MASS 



photometry (ILacy et al. 



photometry, IRAC bands 1 and 2 (3.6 fj,m and 4.5 /xm, respectively) 



20051 ) and our Spitzer IRS mid-IR spectra and MIPS photometry. 



This additional photometry used in the raodelin g procedure is presented in Table [H The 
modehng code is described by lSajina et al.l (120061 ). Modifications to the code to enabl e mod- 
eling of the type-2 objects in our sample were briefly described by iLacy et al.l (120071 ). Here 
we provide a fuller description of the modeling, and the additional modifications required 
to fit the type-1 objects. The modehng is phenomenological in nature due to the current 
lack of a good physical understanding of the nature and geometry of the dust emitting re- 
gions, particularly in the mid-IR, thus many of the parameters do not have a direct physical 
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Fig. 1. — Left IRS spectra for our type-1 quasar sample. Right IRS spectra for our type-2 
quasar sample. The spectra have been smoothed with a gaussian with sigma of five pixels. 
Vertical dotted hues indicate the wavelengths for 6.2, 7.7 and 8.6 //m PAH features, and 9.7 
and 18 /im silicate features. 
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interpretation. Nevertheless, we considered the effort worthwhile as it allows an empirical 
breakdown of the contributions of starlight, hot dust emission from the AGN, and cooler 
dust emission. 

The type-2 SEDs were modeled using the following components. (1) A single stellar 
population from iBruzual fc CharlotI (120031 ) with an age of 5 Gyr to represent the evolved 
stellar population. (2) A power-law component for the mid-IR emission, with an exponential 
cutoff at short wavelengths to represent dust sublimation, and a Fermi function cutoff at 
long wavelengths designed to match the SED shape of the least star-forming object in the 
sample of type-2s, SST 1721-1-6012. The functional form used was: 



-^AGN^ 

e(i'-i'hcut)/«' -|- 1 



where the normalization, I/agn; power-law index a and Tsub (a proxy for the subli- 
mation temperature) are allowed to vary, and i^hcut and w were fixed at 0.11 x lO^^Hz and 
0.017 X lO^^Hz, respectively {h and k-Q a re the Planck and Boltz man constants, respectively). 
The Galactic Center extinction curve of IChiar &: Tielend (120061 ) was then applied to this con- 
tinuum to fit any silicate absorption feature. (3) A warm dust component to represent the 
small grain emission from Hll regions, represented by a power-law with fixed high and low 
frequency cutoffs: 



^SG 



^/'^sglg-i^sgh/l^ 



(2) 



This component is poorly constrained in most of the fits, so 7 was fixed at a typical 
value of two (e.g. Sajina et al. 2006). The high and low frequency cutoffs, Uggi and Usgh were 
set to 0.062 X lO^^Hz and 0.17 x lO^^Hz, respectively. (4) a thermal modified blackbody 
model for the far-infrared emission: 



^FIR 



3+/3 



(3) 



For most objects, we fixed the temperature of this component, T, to 45K, because of 
the uncertai nty in our 160 /im photometry. (3 was fixed at 1.5. In addition, the PAH model 
described in lLacy et al.l (120071 ) was added when evidence of PAH emission was seen, and SST 
1721-1-6012 had a silicate emission profile added, whose shape was approximated by that of 
the extinction curve. 



For the type-1 quasars, component (1) was replaced by an SDSS composite quasar 
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spectrum. As the composite of IVanden Berk et al.l (120011 ) suffers from noticeable host galaxy 



contamination at the long wavelength end, we constructed a new composite by subtracting 
the continuum from the composite of Vanden Berk et al., and adding the residual emission 
li ne composite to a new c ontinuum constructed using line-free SED points from the composite 



of [Richards et al.l (120061 ). 



Most type-1 objects also required an extra near-infrared component to match the SEDs. 
This "very hot" component was modeled as a 1000 K modified blackbody (Eq. [3] with 
T = 1000). The physical origin of this component is unclear, bu t it is most like l y due 



to hot dust close to the sublimation radius (e.g. Netzer et al. 2007). iGallagher et al.l (120071 ) 
show that the strength of this component corre lates with quasar lumin osity. We do not find 



this to be true for our sample, but note that IGallagher et al.l (120071 ) have a much larger 
sample size spread over a range of luminosities. 

The modeling results for the type-1 objects are presented in Fig. [21 and the results for 
the type-2 objects are presented in Fig. [31 Several differences between the type-1 sample and 
the type-2 sample are immediately apparent from the fitting results. The type-1 objects have 
a strong optical quasar component that can extend significantly into the mid-IR. The 1000 
K component is significant in all six type-Is and also contributes to the mid-IR. Three of the 
six type-1 object SEDs showed no evidence of the warm small grain dust power-law emission. 
The type-2 object SEDs are dominated in the mid-IR by hot dust with no contribution from 
the optical stellar light and no evidence for the 1000 K component. Also present in the 
mid-IR portion of the type-2 objects is a variety of PAH emission and silicate absorption. 
Five of the six type-2 objects have significant contributions to the mid and far-IR from the 
warm small grain dust power-law. SST 1721+6012 was the only type-2 not fit well with this 
component. 

We note that the results of the fitting procedure for the object SST 1724+5915 under- 
estimate both the observed amount of PAH emission and the far-IR emission. This is likely 
due to the large amount of sihcate absorption also present in the spectrum of the object. In 
particular the 9.7 /im absorption feature can affect the shape of the 7.7 fim PAH feature, 
making it difficult to fit with an empirical template. We include the data in our subsequent 
analysis, but note that star formation rates derived from the far-IR emission and the PAH 
emission likely underestimate the true star formation of this object (see §4.2p . 

We averaged the modeled type-1 SEDs and type-2 SEDs to highlight the differences 
between the two samples. We first normalized the SED of each target at 24 /im to ensure 
each object had equal weight. By chosing to normalize at 24 /im we can compare the relative 
amounts of hot and cool dust in the averaged SEDs. Though at 24 /im, extinction is still 
present, it is significantly reduced compared to shorter wavelengths. It is clear from Fig. H] 
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Fig. 2. — The type-1 sample. The spectra are presented in their rest frame. From top to bot- 
tom are SDSS 1711+5841, 1711+5855, 1713+5950, 1717+5930, 1717+5948, and 1718+5849. 
The photometry measurements are plotted as circles and the IRS data are plotted as a solid 
line. We show error bars for the MIPS 70 /im and 160 fim photometry. Error bars on other 
photometry measurements are smaller than the plotting symbol. The individual components 
of the model are plotted as dot-dashed lines and the overall fit as the solid black line. The 
blue line represents the quasar composite component, the cyan line the 1000 K component, 
the orange line the mid-IR power-law, the green line the warm small grain dust power-law, 
and the red line the cool 45 K modified blackbody. 
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Fig. 3. — The type-2 sample. The spectra are presented in their rest frame. From top to 
bottom are SST 1711+5904, 1711+5858, 1713+5855, 1718+5953, 1721+6012, 1724+5915. 
The plotted symbols are the same as in Fig. [2l The photometry measurements are plotted 
as circles and the IRS data are plotted as a solid line. We show error bars for the MIPS 70 
/xm and 160 /im photometry. Error bars on other photometry measurements are smaller than 
the plotting symbol. The individual components of the model are plotted as dot-dashed lines 
and the overall fit as the solid black line. The blue line represents the host galaxy stellar 
component, the orange line the mid-IR power-law, the green line the warm small grain dust 
power-law, and the red line the cool 45 K modified blackbody. 
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log(A [fj.m]) 

Fig. 4. — Averaged modeled SEDs of type-1 (blue) and type-2 (red) quasars. The spectra 
are normalized at 24 /im. We plot z/f^^ to highlight the differences between the two. 
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that on average type-2 objects show more PAH emission, more sihcate absorption and more 
far-IR emission than type-1 objects. In addition, while the very hot 3 fim bump is clear in 
the type-1 spectrum, no bump is seen in the type-2 spectrum. We note the caveat that the 
type-1 and type-2 objects could have equivalent amounts of PAH emission if that emission 
were significantly diluted by a strong mid-IR continuum in type-1 objects. 

The important question is if these observed differences could arise solely due to orien- 
tation effects. In that scenario the optical through mid-IR light is absorbed and re-radiated 
in the far-IR. To test this, we integrated the optical through far-IR fluxes of the normalized 
and averaged spectra and found them to be identical. These results are consistent with the 
orientation hypothesis, though the geometry and distribution of the obscuring dust cannot 
be proved from this simple test. 

Furthermore, under the orientation hypothesis, the dust tori of the type-1 and type-2 
quasars should have the same intrinsic luminosities. Our modeling results show that the tori 
of the type-2 quasars are slightly more luminous than the tori of the type-1 quasars, though 
the difference is not statistically significant. This may be a result of our mid-IR selection 
technique, as those wavelengths are still somewhat affected by extinction. In the next section 
we examine the PAH emission and the star formation contribution to the far-IR emission to 
gain more insight into whether the far-IR emission is different between the two objects. 

4.2. Star Formation Rates 

We characterized and compared the star formation rates of the quasar host galaxies 
by comparing the PAH emission with the far-IR luminosity of our sample (Fig. We 
measured the far-IR luminosity from only the modeled warm (Eq. [2]) and cold dust (Eq. [3]) 
contributions to the SEDs (hereafter L'^j-^), excluding dust heated by the central AGN and 
emitting in the mid-IR. This definition of Ij'pjfi specifically excludes dust heated by the AGN 
while still measuring dust heated in star forming regions, even if the two emit at similar or 
overlapping wavelengths. This is in contrast to using a defined wavelength cutoff, which 
would measure dust heated by either source. However, it is still important to note that the 
far-IR emitting dust could also be heated by the AGN if that dust were at large enough 
distances from the nucleus. It is important to consider this potential AGN contribution to 
the far-IR emission when calculating SFRs. 



The PAH emission is an independent indicator of the star f ormation a ctivity in the 
host galaxy. UV and X-ray emission from the AGN destroy PAHs (IVoitlll992l ). but soft UV 
emission from star formation excites PAHs into emission. In Fig. [S] we find a correlation 
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Table 1. Targets and their Fluxes in mJy 
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Table 2. Additional Photometry Used in Modeling Procedure 
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19.36 


18.85 


18.92 


18.72 


18.82 


18.53 




16.19 


0.46 


0.61 


SDSS 171748.43+594820.6 


18.42 


18.06 


18.07 


18.14 


18.09 


18.04 




16.05 


0.65 


0.90 


SDSS 171818.14+584905.2 


19.17 


18.70 


18.70 


18.57 


18.70 








0.68 


0.87 


SST 171106.8+590436 


22.30 


21.52 


20.07 


19.51 


19.15 


18.0 


17.2 


16.3 


0.32 


0.44 


SST 171147.4+585839 


22.94 


22.30 


21.61 


20.80 


20.31 


18.8 


18.0 


16.6 


0.52 


0.82 


SST 171324.1+585549 


22.82 


21.83 


20.95 


20.21 


20.06 


19.0 


18.2 


17.3 


0.20 


0.43 


SST 171831.7+595317 


22.93 


21.91 


20.93 


20.18 


19.78 


18.7 


18.0 


17.2 


0.16 


0.26 


SST 172123.1+601214 


21.69 


20.28 


18.99 


18.74 


18.33 


17.8 


16.5 


15.2 


1.03 


1.60 


SST 172458.3+591545 


21.94 


21.69 


20.31 


19.61 


19.17 


18.1 


17.2 


16.3 


0.43 


0.49 



Note. — Additional photometry for the modeUng procedure was retreived f rom NASA Extragalactic Database 
(NED). SDSS u, g, r, i, and z magnitudes are listed in the AB system llOkelll990ll . The u and g photometry were not 
used in the fitting of the type-2 targets. 2MASS J, H and Ks magnitudes are presented in Vega system. IRAC bands 
1 and 2 (3.6 lira and 4.5 fan) are listed in mJy. 
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12.0 
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Fig. 5. — Top: Integrated PAH luminosity as measured from our scaled template is plotted 
against far-IR luminosity. Bottom: Relative amount of PAH emission compared to the far-IR 
dust emission. In both plots, L'^^^ was measured from the warm and cold far-IR components 
of the model only. Lpah and L'pj^ are in Lq units. Type-1 quasars are plotted as diamonds, 
and type-2 quasar s are plotted as Xs. T he starburst galaxies NGC 7714 and NGC 2623 are 
plotted as circles (IMarshall et al.l 120071 ). The dashed line indicates the median ratio of the 
type-1 objects of our sample, and the dotted line indicates the median ratio of the type-2 
objects. 
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between PAH emission and far-IR emission similar to previous works (e.g., Schweitzer et 
al. 2006, Menendez-Delmestre et al. 2008). We also normalized the PAH emission to the 
far-IR emission, so we can compare the relative contribution of star formation to the far-IR 
emission. We measured the PAH luminosities by integrating our scaled templates, including 
the full profiles of each individual PAH feature in the template. This is different from many 
previous works (e.g., Schweitzer et al. 2006, Menendez-Delmestre et al. 2009), where only 
individual PAH features are measured, and definite wavelength cutoffs are applied. We find 
that the quasars as a whole have a fairly constant ratio across a range of far-IR luminosities. 
The median Ijpah/^'fir ratio for the type-1 objects of our sample is 0.012, and the median 
for the type-2 objects is 0.021. We note that the difference between the type-1 objects 
and type-2 objects is not significant consid ering the scatter in th e sample. The median of 



our overall sample is 0.017. In comparison, [Marshall et al.l (120071 ) measure PAH luminosity 



using a similar method to ours and find that the starburst galaxies NGC 7714 and NGC 2623 
have 'LpAH I^FiR ratios of 0.052 and 0.021, respectively. We cannot draw strong conclusions 
based on measurements of two objects, and note that many objects of our sample have ratios 
consistent with these objects. However, the average Ijpah/^fir from Marshall et al. is 0.037, 
higher than the median value of our sample. If taken at face value, these numbers indicate 
that some of the cool far-IR emitting dust is heated by AGN emission instead of stellar 
emission. 

In Table |3] we present the integrated mid- to far-IR emission (8 — 1000 /im, hereafter 
L/_r). We find that nearly all of our targets have log(L/R/L0) > 11.75. SST 1721+6012 
is the only object with XogihiR/ho) < 11.75. We also hst the star formation rates (SFR) 
of our sample in Table [31 We calculated SFRs in two different ways. First, we converted 



the in tegrated warm and cold dust components of our model into SFRs using the iKennicutt 



(119981 ) calibration for starburst galaxies: 



SFR [Mq yr-^] = 4.5 x lO'^^L^^^ [erg s'^] (4) 

Although our objects have IR luminosities similar to those of ULIRGs, their SFRs are 
lower than those of typical ULIRGs. This is because some of the IR luminosity of our objects 
is due to AGN-heated dust, and we explicitly exclude this component when calculating SFRs. 
The SFRs of our sample are instead typical of LIRGs, which can range from 20 — 200 Mq yr~^ 
(e.g., Sanders et al. 2003). We find that the type-2 objects show higher SFRs on average (92 
Mq yr~^) than the type-1 objects (67 Mq yr~^). However, the scatter is large, and the object 
with the lowest SFR is a type-2, SST 1721+6012. We also independently calculated the SFR 
from the PAH emission. This can be done by predicting a far-IR luminosity based on the 
PAH luminosity. We integrated our scaled PAH template to measure the PAH luminosity 
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of each objec t in our sample . We adopt an average ratio of Lpah/^'Fir ~ 0.037 based on 
the results of [Marshall et al.l (120071 ). Subsequently, for all objects in our sample except two, 
the SFR predicted from the PAH luminosity is lower than the SFR predicted by integrating 
the warm and cold dust components of the far-IR (L^^^) . Again we found that the average 
type-2 quasar has a higher SFR (87 M© yr~^) than the average type-1 quasar (37 Mq yr~^). 



5. Summary and Discussion 

We obtained mid-IR spectra and far-IR photometry for six type-1 and six type-2 mid-IR 
selected quasars. We modeled the SEDs of our sample from the optical to the far-IR using 
available SDSS, 2MASS and IRAC photometry, in addition to our IRS spectra and MIPS 
photometry. We compared the modeled SEDs of the quasars and calculated SFRs based on 
PAH and far-IR emission. Our main results are: 

(1) The type-1 quasar mid-IR spectra show a featureless continuum with significant 3 
Hm emission. The full optical through far-IR SEDs are typically fit well with an optical 
quasar component, a 1000 K modified blackbody, a mid-IR power-law (~ 3 — 30yum), and 
a cool dust far-IR modified blackbody (45 K). Three of the objects also show evidence of 
warm small grain dust fit with a power-law (~ 25 — 65/im). 

(2) The type-2 quasar mid-IR spectra show a variety of PAH emission and silicate 
absorption strengths. The full optical through far-IR SEDs are fit well using a stellar pop- 
ulation in the optical, a mid-IR power-law (~ 3 — 30/im), varying amounts of a PAH emis- 
sion template and silicate absorption in the mid-IR, warm small grain power-law emission 
(~ 25 — 65/im) as well as a cool dust modified blackbody (45 K). 

(3) Our averaged and normalized modeled spectra show that type-2 quasars exhibit 
more far-IR emission than type-1 quasars. The excess is roughly equivalent to the deficit 
in the optical with respect to the type-1 quasars, which is consistent with the orientation 
hypothesis. In this scenario the increased far-IR emission arises from obscuring dust that is 
reradiating absorbed light. 

(4) We measured Lpah by integrating our scaled template for each object. We compared 
the typical Lp^H to L^^^ ratio of our sample to that of starburst galaxies that have had Lp^H 
measured in a similar manner as the method we used. Within our sample, we found that 
the type-2 objects have a larger ratio than the type-1 objects. We note that the difference 
between the two is not significant. Con sidering our overall sa mple, we found that the ratio is 



smaller than for the starburst galaxies (IMarshall et al.ll2007l ). This result suggests some far- 



IR emitting dust can be warmed by the central AGN rather than by young stars. However, 
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Table 3. Luminosities and Star Formation Rates 



Target 


log(L,^/Lo)" 


log(LpAH/L0)^ 


SFR [Mo yi-^Y 


SFR [Mq 


SDSS 171117.66+584123.8 


12 


25 


9.3 


109 


10 


SDSS 171126.94+585544.2 


11 


82 


9.4 


40 


12 


SDSS 171334.03+595028.3 


11 


95 


9.3 


27 


10 


SDSS 171736.91+593011.5 


12 


05 


9.3 


37 


9 


SDSS 171748.43+594820.6 


12 


06 


10.3 


66 


98 


SDSS 171818.14+584905.2 


12 


26 


10.2 


125 


82 


SST 171106.8+590436 


11 


79 


10.0 


73 


48 


SST 171147.4+585839 


12 


46 


10.9 


183 


335 


SST 171324.1+585549 


12 


13 


10.0 


85 


52 


SST 171831.7+595317 


12 


47 


10.0 


143 


50 


SST 172123.1+601214 


11 


48 


9.0 


14 


4 


SST 172458.3+591545* 


11 


77 


9.8 


54 


29 



'^The luminostiy L//j corresponds to the integrated modeled SED from 8-1000 /xm. 

'^The PAH luminosity is measured from the integrated profiles of all PAH features in the 
model. 

'^Star formation rates derived from the sum o f the warni and c old dust components of the 



model only. The SFR was calculated using the iKennicuttI (119981 ) relation. 



'^Sta r formation rates derived from PAH l uminosity using t he results of [Marshall et al. 



( 120071 ) to convert from hp ah to hpiR and the IKennicuttI (119981 ) relation 



*The fit of SST 1724+5915 underestimates the PAH luminosity and far-IR emission, thus 
also underestimating the SFRs derived from these quantities. 
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a more direct comparison between AGN and starburst galaxies using larger samples would 
be needed to further quantify the amount of AGN-heated far-IR emitting dust. 

(5) We calculated the SFRs of our sample using two methods. First we integrated the 
warm small grain and cool modif ied blackbody era ission (excluding all other components) of 



converted the Lpah to Lpm (jMarshall et al. 



the modeled SEDs and used the iKennicuttI (1998 ) relation to calculate the SFRs. We also 



20071 ) to measure the SFRs based on the PAH 



emission only. In each case the type-2 objects had higher SFR on average than the type-1 
objects. However, the difference is not significant given the scatter. 

(6) We integrated our modeled SEDs from 8 — 1000 /im and found that seven of the 12 
can be considered ULIRGs (log(L/i^/L0) > 12), while four others have log(L/ij/L0) > 11.75. 
The sole exception is SST 1721+6012, which has log(Lj-/j/L0) = 11.48. While our sample 
appears to have IR luminosities similar to ULIRGs, the quasars exhibit SFRs more typical 
of LIRGs. We found the AGN can heat most of the dust in the mid-IR and even some in 
the far-IR. Thus measuring SFR in quasars from the full 8 — 1000 /im SED overestimates 
the amount of dust warmed by young stars and the true SFR. 

(7) The modeling results suggest that type-1 quasars can have mid-IR fluxes that have 
significant contributions from the optical power-law and 3 /im bump. This 3 /im component is 
not required for the type-2 SEDs, and the hot dust most likely originates near the sublimation 
radius of the AGN. If this component is present in the type-2 quasars, it must be at least 



partia lly obscured by cooler dust, either a dusty torus or dust in the host galaxy. iLacy et al. 



(120071 ) showed that the host galaxies of the type-2 quasars have a range of inclinations. 
Therefore, the obscuration of the 3 /im bump, if it is intrinsically present, should be from 
material close to the nucleus and independent of galaxy orientation. 

While our analysis does not distinguish definitively between orientation and evolution 
based unification schemes, it does directly address the dust content, ongoing star formation in 
the host galaxies and the AGN contribution to the far-IR emission. Our results are consistent 
with the orientation hypothesis, but we cannot rule out an evolutionary connection. The star 
formation based on PAH measurements is a larger contributor to the far-IR luminosity than 
the AGN for type-2 objects compared to type-1 objects. This difference is only suggestive 
due to the scatter of both Lpah and L^^^ in our sample and could arise from intrinsic 
differences in the bolometric luminosities of the objects. The Lpah/L'pjj^ ratio in both type- 
1 and type-2 quasars is lower than that of starburst galaxies, suggesting that some emission 
in the far-IR is from dust heated by the AGN and not star formation. This is an area for 
continued research, as our analysis was based on few objects. 
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